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Abstract. The continuum Kardar-Parisi-Zhang equation in one dimension is lattice 
discretized in such a way that the drift part is divergence free. This allows to 
determine explicitly the stationary measures. We map the lattice KPZ equation to 
a bosonic field theory which has a cubic anti-hermitian nonlinearity. Thereby it is 
established that the stationary two-point function spreads superdiffusively. 



1 Introduction 



Oversimplified models of surface growth have remained of interest, in particular 
their one- dimensional version. In very general terms, one investigates a conserved 
field with a stochastic dynamics which does not satisfy the condition of detailed 
balance. Mostly models with discrete space and discrete heights have been in fo- 
cus, as e.g. the asymmetric exclusion processes (ASEP), also its totally asymmetric 
version (TASEP), and the polynuclear growth (PNG) model. We refer to [H |2l |3] 
for reviews. Some features of these models are related to integrable systems and 
random matrix theory, which has become an independent motivation for intense 
study. There are unexpected connections to random tilings, as the Aztec diamond, 
and to crystal shapes, somewhat less surprising to directed last passage percola- 
tion, also known as directed polymer in a random medium. For several features of 
ASEP, TASEP, and PNG exact solutions can be obtained through a rather intricate 
asymptotic analysis. In particular, the stationary two-point function has been com- 
puted in the scaling limit [H [5] . Except for model dependent scale factors one finds 
the same scaling function in the various models, supporting, at least partially, the 
universahty hypothesis for growth models with local deposition rules and neglected 
surface diffusion. 

In a seminal paper, Kardar, Parisi, and Zhang [6] proposed a continuum equation 
for surface growth. In one dimension their equation reads 



Here h is the height function over R at time t, Ab is the strength of the nonlinear 
growth velocity, z/b is the strength of the local smoothening, and ^ is normalized 
Gaussian space-time white noise, {^{x,t)^{x',t')) = S{x — x')5{t — t'). The slope 
u = dxh then satisfies the stochastic Burgers equation. 



Note that this relation is a special property of one dimension. 

It is expected that Eq. (11. 2p is in the same universality class as ASEP and 
PNG. The most convincing evidence comes from the numerical simulation of the 
discretized version of (II. ip with 1024 lattice sites The stationary two-point 
function is computed and good agreement is obtained with the exact solution from 
TASEP and PNG over a substantial range of wave vectors. 

On the theoretical side, one notes that the noise in (II. 2p is very singular. To 
cope with this difficulty, one worked out approach [8J is to make Eq. (II. ip linear 
through the Cole-Hopf transformation 



dth{x, t) = |Ab(9,/i(a;, t)f + u^dlh{x, t) + ^/l\i{x, t) . 








Then the "partition function" Z{x,t) satisfies 



hZ{x,t) = u^dlZ{x,t) + 



2z/b 



^{x,t)Z{x,t) 
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with initial conditions Z{x,0) > 0. Being linear, one can give sense to Z{x,t) as 
a stochastic process with continuous sample paths. In particular, Z{x,t) > with 
probability 1 and one defines h{x, t) through Eq. (11.31) . In fact, this type of solution 
of the KPZ equation can be recovered by a suitable continuum limit of the ASEP 
with a properly chosen weak asymmetry. 

An alternative approach is to regularize the noisy Burgers equation (II. 2p . De- 
pending on the point view there are then two choices. The fluid dynamics camp 
regards (1.2) as oversimplified large scale randomly stirred Navier-Stokes equations 
and studies small scale properties of the solution, see [TOl [11] out of a large body 
of literature. But then it is natural to keep the 5-function in time and to replace 
the 5-function in space by a smoothened version. On the other hand the surface 
growth community knows that the lattice structure of the solid defines a smallest 
length scale. Thus one discretizes the KPZ equation, as done without further ado 
in any numerical simulation, and studies its properties on large scales. In relation 
to ASEP and PNG it is natural to follow the second route and to discretize M as 
5Z with lattice constant 6 > 0. The field variables are denoted by Uj{t) G M, j G Z, 
as the lattice version of u{x,t) in (11.21) . Then the lattice KPZ equation, studied in 
this paper, reads 



Here € Z} is a collection of independent and normalized white noises. Below 

we will determine a family of stationary measures of (II. 5p . We consider then the 
space-time stationary solutions to (II. 5p and the object of prime interest will be 
stationary two-point function 



at average slope p = {uj{t))p. 

Each one of the three terms on the right hand side of (II. 50 converges as 5 — to 
the corresponding term in (II. 2p . To obtain the scale invariant theory, however, one 
has to adjust the "coupling constants" in (II. 5p in such a way that S{j,t) converges 
to a nontrivial limit. This is the reason for the subscript b, which stands for bare 
coupling parameters. 

To provide a brief summary, in Section [2] we rewrite (II. 5p as a bosonic field 
theory with a cubic anti-hermitian nonlinearity, where for simplicity we consider 
only the case of zero slope, p = 0. As explained in Section [3], restricting the bosonic 
occupation variables to 0,1, i.e. restricting to hard core bosons, yields precisely the 
stationary ASEP at density 1/2. The bosonic representation is used in the study of 
the spreading of Sp=o{j, t). Sp=o{j, t) is even in j and the normalized variance reads 




|(Ab/35)(M^+i + - - 

+ (z/b/5')(M,+i - 2u, + n,_i) + (Db/5)'/' 



J el. (1.5) 



(1.6) 




(1.7) 
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In Section [5] we will establish the upper and lower bounds as 

t^/^ < Var(t) < t^/^ 



(1.8) 



as t — s> oo. We then develop an iterative scheme based on the relaxation time 
approximation and recover the KPZ prediction of Var(t) = f^^^. Very recently J. 
Quastel announced a proof of the following bounds, 

c_t^/3 < VarcH(t) < c+t^^^ (1.9) 

for large t and suitable constants < c_ < c+ [9]. Here VarcH(^) is computed 
from the variance of log Z{x,t) with Z{x,t) the Cole-Hopf solution of ( 11. 4p with 
initial data Z{x,0) such that logZ(x, 0) is two-sided Brownian motion in x. In the 
final section the continuum approximation is discussed and contrasted with the well 
studied case of (11. 2p for D^, = 0, z/b ^ 0, and white noise initial data. 

2 The lattice KPZ equation 

To simplify notation we set Aq = Ab/5, i^o = ^h/^"^, Dq = D\,/5, and a = uq/ Dq. 

The lattice KPZ equation (11.51) conserves the slope field u, which becomes man- 
ifest through introducing the current function 

= l^oiu^j + UjUj+i + M^+i) + iyo{uj+i - Uj) 

= Wj + uo{uj+i - Uj) . (2.1) 

Then ^ 

—Uj = Wj - Wj-i + a/A)(^j - ^j„i) , j eZ. (2.2) 

In principle, there are many possibilities to discretize d^ulx)"^ in (II. 2p . Our choice 
is singled out by the facts (i) the discretization involves only nearest neighbors and 
(ii) the drift term is a divergence free vector field, i.e. 

J2djiwj-Wj.^) = 0, (2.3) 

where dj = d/duj. For this particular discretization one can compute explicitly the 
invariant measures. (12. ip and (12. 2p was first proposed in p^, see also |13] for a more 
recent study. 

For a while we study (12. 2p on the ring [1, . . . , A^], i.e. 

UN+j = Uj and ^n+j = • (2-4) 
The generator for the diffusion process (12.21) is then 

N 
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Noting the identity 

N 

^{wj - Wj-i)uj = , (2.6) 
i=i 

one verifies that, for every p G M, the independent Gaussians 

N 

n{(«/vr)^/^exp[-a(^, -p)^]} = (2.7) 

are invariant for (12.51) . i.e. for every smooth function / on configuration space it 
holds 



Jm.n 



Note that the stationary measure does not depend on Aq. Let us denote the average 
with respect to (^^p)^ by {■)p,N, {^)p,N = 1- Then 

{uj)p,N = P (2.9) 

is the average slope, 

(«?)p,7v - {u,)In = ^ (2.10) 

the variance, and 

{w,),,N = lXo{a~' + 3p')=jip) (2.11) 

the average current. For simplicity we will consider only the slope case, p = 0, 
and hence omit the index p. 

We switch from the generator Ln to the hamiltonian Hn through the ground 
state transformation 

^P^LNi^^)-' = -HN. (2.12) 

Then 

N N 

Hn = lDoJ2i - i^j - 5j-i)' + a\uj - Uj^iY - 2a) - J^i^j - Wj^i)d, . (2.13) 
i=i i=i 

We introduce the standard annihilation/creation operators at site j through 

1 1 

{auj + dj) , a* = {auj — dj) . (2-14) 



^ V2a ^ V2a 



a 

Note that they satisfy the canonical commutation relations 

[tti, a*] = 6ij . (2.15) 

Then 

= Ho^N + Xo{A% - An) , (2.16) 
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where 

N 

Ho,N = T^o ~ - %) (2-17) 

and 

N 

ijv = (3 • 2^/^)-^q;-^/2 ^ (aja*+iaj+i + %%a*+i - 0*0^+10^+1 - 0*0^0^+1) . (2.18) 

We want to compute the propagator exp[— iifjv], t > 0. To reduce the number 
of parameters we rescale time as t -w jjQt. Then the prefactor of Hq j^ becomes one. 
We also introduce the couphng constant 

A = (3-2=^/2)-%q;-VV- (2.19) 
Hn = Ho,N + KA*j, - An) , (2.20) 

N 

Ho,N = ~ %)*(%+! - %) ' (2-21) 

AT 

Ajv = ^ + ~ ~ Qj^^J+i) ■ (2-22) 

i=i 

It is convenient to use the Fock space representation of the bosonic field aj,a*. 
Physically this corresponds to a study of local excitations away from the Gaussian 
stationary measure, as e.g. encoded by the two-point function S{j,t). For N — 00 
the Fock space ^ is over £2 — •^2(Z) with the n-particle space — (•^2)^m 



Then 
where 



^ = 0S^n. (2.23) 



n=0 

An element / e ^ is a sequence {/o, /i, ...,/„,.. .}. /„ : Z" — > C and /„ is symmetric 
in its arguments, /o G C. The scalar product on is 

{f:9)n= ^ fn{Xl,...,Xnygn{Xu...,Xn). (2.24) 

/ e 5^ if and only if 

oo 

(/,/) = I/or +E(/-/")"<°°- (2-25) 

n=l 

is the n-particle subspace. The Fock vacuum is the vector fl — (1, 0, 0, . . .). 
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The {aj,a*,j G Z} are the usual bosonic annihilation/creation operators on ^J, 
which are defined through 

{ajfn+i){xi, ...,Xn) = Vn + ...,Xn) (2.26) 

and 



{a*fn-l){Xi, ...,Xn) 



1 " 

-^y^S{xe- j)/„_i(xi, ...,xi,...,Xn) (2.27) 



with the convention that Xi means the omission of X£ from the configuration (xi, 
. . . , Xn), see e.g. ^A\. From fl2.20p - (12.221) in the hmit oo, we arrive at 

H = Ho + X{A* - A) . (2.28) 

Exphcitly 

Ho = ^(flj+i - aj)*(aj+i - aj) (2.29) 

and 

A = (aja*^;^aj_|_i + ajaja*_,_]^ — a*aj+iaj+i — a*ajaj+i) (2.30) 
jez 

as operators on Fock space. 

Ho \ is the Laplacian and corresponds to n independent random walks on 
Z with nearest neighbor hopping of rate 1. Clearly, Hq = Hq, while A* — A is 
antisymmetric, {A* - A)* = -{A* - A). Using ([226]), (HYlh . the action of A on 
Fock space is represented by 



Afn+l{Xi, ...,Xn) = + 1 (/„+i(xi, . . . ,Xn,Xe - 1] 



£=1 

~fn+l{Xi, . . . , Xn, Xi + 1) + fn+l(yXi, . . . ,Xi, . . . , X„, Xi — l,Xi — 1) 

-/«+i(a;i, + l,a;^ + 1)) . (2.31) 

To have an example, let us consider the two-point function at slope zero, 

Sij,t)=E{u,it)uom- (2.32) 

Here N = oo and the expection E(-) refers to the solution of (12. ip . (12. 2p with Uj{0) 
distributed as independent Gaussians with mean and variance l/2a. In Fock space 
S{j, t) is expressed through 

5(j, t) = — (n, aoe-^o*^a*n) , t>0. (2.33) 
2a ■' 

Note that by the refiection Uj ^ —U-j and by time stationarity, one concludes 

SU,t) = S{-j,t), S{j,t) = SU,-t), t>0. (2.34) 
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Computationally it is often more convenient to switch to Fourier space. For 
/ : Z ^ C we define the Fourier transform 

/(A;) = 5^e-^2-^V(j) (2.35) 

with inverse ^ 

fU) = [ dfce'2-^V(fc) • (2.36) 

Thus the one-particle space is now L'^{[0, 1], dk). We set T = [0, 1] as first Brillouin 
zone. Correspondingly 

a{k) = J2^"^"^^(^j- (2.37) 

Note that 

[a{k),a{k')*] = 5{k-k'), (2.38) 

which means that the corresponding map between Fock spaces is unitary. In Fourier 
representation it holds 



Ho = / dkuj{k)a{k)*a{k) , u{k) = 2(1 - cos(27rA;)) (2.39) 
Jj 

and 

A = —i dkidk2dksS{ki — k2 — k^)(2sm{27Tki) 

+ sin(27rA;2) + sin(27rA;3))a(A;i)*a(A;2)a(A;3) . (2.40) 



3 Relation to the ASEP at density 1/2 

The partially asjTumetric simple exclusion process (ASEP) is a stochastic particle 
system on Z with hard exclusion, i.e. the occupation variable rjj at site j takes only 
the values 0,1. The state space is {0, 1}^. Particles hop to the right neighbor with 
rate 1+p and to the left neighbor with rate 1 — p, \p\ < 1, provided the destination 
site is empty. Therefore the generator. Las, reads 

LAsfiv) = E ((1 + - ^^+1) + (1 - - vj)m+i) ifiv''^') - fiv)) , (3.1) 

where rj^^'^^ denotes the configuration rj with the occupations at sites j and j + 1 
interchanged. The Bernoulli measures are invariant under Las. We consider only 
the stationary process with Bernoulli 1/2. 

We expand Las in the natural basis for the Bernoulli measure with density 1/2 
|15j . The basis functions are labelled by finite subsets. A, of Z and are of the form 

MV) = n(2^^- - 1) • (3-2) 
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In particular, with (■)i/2 denoting average over Bernoulli 1/2, 

(^AiV^A2)l/2 = 5(Al,A2). (3.3) 

A general function, v, is represented by 

t;(r/) = 5^t;(A)^A(r7), ^ |{)(A)|2 < oo , (3.4) 

AcZ AcZ 

the sum being over all finite subsets of Z. In this basis Las is represented by 

Las = -Has, Has = S + piA* - A) . (3.5) 

Here S = S* and 

SviA) = -J2 (^(A.,.+i) - v{A)) , (3.6) 

where A^^^^i is obtained from A by exchanging the occupancies at x and x + 1. 
To define A and A* we introduce the outer left boundary of A, £(A) = {x\x ^ 
A,x + 1 G A}, and the outer right boundary of A,r(A) = {x\x ^ A,x — 1 G A}. 
Correspondingly, the inner left and inner right boundary of A are i{A) = {x\x G 
A, X - 1 ^ A}, f(A) = {x\x G A, X + 1 ^ A}. Then 

Av{A) = ^ v{AU{x})- ^ v{AU{x}), 
xeeiA) x-er(A) 

A*v{A) = J2v{A\ {x}) - {^}) • (3-7) 

xee{A) xef{A) 

Las is the generator in the number space representation. To be able to compare 
with H one still has to transform unitarily to Fock space representation. For this 
purpose let |A| = n, A = {xi, . . . ,Xn}- Out of '0(A) we construct G dn by 

/'n('^l) • • • ) ^n) 



{1/ Vn\)v{{xi, . . . , Xn}) (xi, . . . , Xn) has no coinciding points , 
otherwise . 

(3.8) 

We observe that 



(/n; fn)n ^ ^ \fn{xi, . . . , Xn) \ 

= 1 5^ \v{{x,,...,x4)\'= Yl 1^(^)1'- (3-9) 

' xeZ",xe^x,n, all ^^m AcZ,|A|=ri 

Hence the map defined by (13. 8p is unitary. The set of all v^s such that the right 
hand side of (13. 9p is finite defines the subspace P^n of dn- Pdn consists of functions 
which vanish whenever their argument {xi, . . . ,Xn) has coinciding points. P is a 
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projection operator on ^„ and P as an operator on ^ is defined through its action 
on each subspace "Sn- 

Now, if fn+i € P^n+i, then in the expression (12.311) for A the two last summands 
vanish. A carries the prefactor + 1 which is balanced by the normalization l/y/n\ 
in (13. 8p . Thus we conclude that, for X = p, 

PHP = Hj,s (3.10) 

on Fock space. Restricting the lattice KPZ hamiltonian to the subspace P^ results 
in the hamiltonian of the asymmetric simple exclusion process. Note that Las is 
the generator of a Markov jump process only for |p| < 1. On the other hand, H is 
defined for all A and thus (13.101) provides an "analytic continuation" of Las- 

It would be most useful to have comparison inequalities between lattice KPZ 
and ASEP. For the symmetric part of the lattice KPZ hamiltonian, the quadratic 
form on ^„ reads 

n 

W- (3.11) 

£=1 xi,...,x,i(^Z 

Thus Hq restricted to is the Laplacian on Z". On the projected operator 
PHqP corresponds to the Neumann restriction of the Laplacian (13. lip to the set 
{x G 7^ Xm for all pairs i ^ m}. Hence for the symmetric part it holds 

Ho > PHoP (3.12) 

on ^. 

On the other hand, there seems to be no such simple relation for the asymmetric 
part. In |T5] the authors approximate the ASEP by a model without hard exclusion, 
in such a way that its restriction to the range of P agrees with the ASEP. Restricted 
to the subspace 0^=0 they bound the resolvent of one hamiltonian in terms 
of the resolvent of the other hamiltonian, and vice versa. Their bounds are not 
uniform in n. Since the approximation in [15] has some similarity to the lattice 
KPZ hamiltonian, we expect that their bounds would carry over to H from (12.280 . 

4 Variance of the two-point function for large t 

Since, by the conservation law, 

X = J2Sij,0) = Y,S{j,t), x = ^, (4.1) 

and since, by (I2.34p . 

J2jS{j,t) = 0, (4.2) 
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it is natural to introduce the normalized variance 

Var(t) = x-^$^f5(j,t). (4.3) 

We first represent S{j,t) through the transition probability e^*, where L is the 
generator Ljv from (12.51) with the sum over all j G Z. Then 

SU, t) = E(Mo(0)M,(t)) = (Moe%) . (4.4) 

The expectation E has been defined below (I2.32p and (■) refers to the average over 
the initial data, i.e. over independent Gaussians with mean and variance l/2a. Uq 
and e^*Mj are functions on configuration space. The conservation law is used again 
to partially integrate twice, 

xVar(t) = ^j2(Moe%) 

= / ds^j2(uoe^^(w;, -w,_i)) 

= / ds5^(-2j-l)((now;,)+ rds'((L*no)e^^'w,)) 
Jo ^ Jo ^ 

= 2tMul)+2 fds rds'^(((w;o- W)e^^'(«;,-K)))). (4.5) 
Jo Jo 

It is convenient to switch to Laplace transform. Then 

POO 

xVa.T{C) = / dte-'^'^°*xVar(t) 



= 2C'i^,'{ul) + 2(z/oC)-' - <^o))^773t(^^- - K))) (4-6) 

for C > 0. 

On the other hand, the KPZ scaling theory asserts that 

S{j,t) - xi2\lxt')-'^'fKPz{i2\lxtT'^'j) (4.7) 

with X defined in (14. ip . the renormalized coupling constant j"{0) = Aq, and the 
universal stationary KPZ scaling function /kpz- The validity of (14.71) has been 
proved for the stationary PNG model [4| and for the stationary TASEP [5j. We 
refer to these papers for the definition of /kpz- In particular, 

Var(t) = {2Xlxt'f'{x')KPz (4.8) 



11 



for t ^ cxD with (x^)kpz = / dxx^/KPz(a^) = 0.510523... a model independent 
number. Hence 

Var(C) = (2A^x)^/^r(7/3)(z/oC)-'/=^(x^)KPZ (4.9) 

for — >■ 0. Comparing (14.91) and (14. 6 p yields the prediction for the small ( behavior 
of the resolvent in (14.61) . But before we have to reexpress this resolvent in Fock 
space. 

Let us define the total momentum operator, Ptot, as 



Ptot = 2n dkka{k)*a{k) . (4.10) 
Jt 

Since H is translation invariant, 

[H,P,,,] = 0, (4.11) 
there exists the direct integral decompositions corresponding to Ptot as 

d= r m{k) (4.12) 



Jt 

and 

H= dkH{k). (4.13) 
Jt 

We will need only the fiber at A; = 0. To construct ^5^(0) we consider / G 5^, in the 
momentum representation, such that fo = and are continuous on T". Then 
/° G diO) is defined by 

/°(A;i, . . . , fc„) = fniku • • • , K)6ik, + ... + kn), n = 1, 2, . . . . (4.14) 

The scalar product is given by 

(/n,/n)n= / dk^ . . . dkj {k, + . . . + k^)\ ^k,, . . . , k^)\\ (4.15) 

Jt" 

oo 

(/°,/T = E(/-/nT- (4.16) 

n=l 

is the completion with respect to this norm. 
The free part on is multiplication by 

n 

Qn{ki,...,k^) = Y,^ike). (4.17) 

e=i 

Hence 

{Hof^){ku . . . , fc„) = n^iki, . . . , kn)fniki, . . . , fc„)5(fci + ... + kn) (4.18) 
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and i/o/° £ 1^(0) • Correspondingly, using ( I2.40p . the cubic operator A on acts as 

{Afn+l)ih,...,kn) 

n „ 

= —2iy/n + 1 f / dA;„+i ( sin(27rA;£) + sin(27rA;„+i)) 

x/„+i(/ci, ...,ki,..., kn+i, ki - kn+i) (4.19) 

with a similar expression for A*, see below. Clearly, if f^+i contains the delta 
function 6{ki + . . . + /c„+i), then Afn+i is proportional to 6{ki + . . . + kn) and hence 
e dn{0). Note that H \ di{0) = and H* f ^^(O) = 0. 

The current function at j = with its average subtracted equals wq — (wq) = 
J^o{ui — uq) + (Ao/6)(-Uq + uqUi + uf — a~^). To construct the corresponding vector 

G ^, we use 

u,^o = (2a)-^/'a*^o , (^^? - (2a)-^)^o = (2a)-^af 7/>o , (4.20) 

and the fact that the ground state ipo is mapped to the Fock vacuum. By using 
fl2.27p one obtains = for n = 0, n > 3, and 

/r(xi) = M2a)-'/^{6{xi - 1) - 5(xi)) , 

f^{x,,X2) = {Xo/6){2V2a)-\26{x,)6{x2) + 26{x, - l)5(x2 - 1) 

+(5(xi)5(x2 - 1) + 5ixi - 1)5(X2)) , (4.21) 

which in Fourier space reads 

frik,) = M2a)-'/\e-'''^'^ - 1) , 

f^(k^^ k^) = (Ao/6)(2V2«)"i(2 + 2e-^2-(^i+'=2) + Q-^^^k^ ^ ^-i2nk,^ _ ^4 ^2) 
ywo g jj^g-j jg i^j^gj^ given by 

r° = Xo{2V2a)-'g' , ^° = (0, 0, w', 0,0,...), 

w°(A;i,fc2) =u'(A;i)5(A;i + A;2), w(A;i) = |(2 + cos(27rA;i)) , (4.23) 

enforcing the normalization w{0) = 1. 

With this input the resolvent from (14. 6 p reads 

E((^o - H))— ^(i^, - {w,))) = A^(2^/2a)-Vo-^^^ (C + if)^Y)° • (4.24) 

Comparing (14. 9 p and (14. 6 p together with (14.240 one arrives at the following prediction 
from the KPZ scaling theory, 

{g', (C + H)-'gy = 3-2/3r(7/3)(x2)KPz(A^C)-'/' (4.25) 

for small (. 
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5 Bounds and relaxation time approximation 

While at present we have no techniques to estabhsh (14.251) . following the methods in 
[15] one can study upper and lower bounds for the matrix element {g^, {( + H)^^g^)^. 
Restricting H to 0^=iS^m(O) yields bounds depending on n, which alternate as 
upper and lower bounds and which converge to a limit as n ^ oo. More precisely, 
we define recursively 

T, = iC + Ho)-\ 

T„ = (C + ifo + X^ATn-iA*)-' , (5.1) 

and set 

6„(C) = (^°,T„^T- (5.2) 

E.g., 

64(C) = {9' ,{C + Ho + X'A{C + Ho + X'AiC + Ho)-'A*)-'A*}-'gy . (5.3) 
Then 



and 



hiO<h{0<...<b,{0<h{C), (5.4) 



(^°,(C + if)-V)°= lim6„(C). (5.5) 



The vector g^ has ^2 = the only nonvanishing entry. Hq preserves the 

particle number, while A decreases it by one and A* increases it by one. Hence it 
will be convenient to introduce a notation displaying these special features. Let P„ 
be the projection onto "^n- We set 

PnHoPn = Ho,n (5.6) 

and 

PnAPn+l = An^n+1 , Pn+lA* Pn = An+l^n ■ (5.7) 

For a given n, we define recursively 

f/t\ = X^Ar,.i,n{C + i/o,n)-'^n,n-l , (5.8) 
f/i") = X^Am,m+l{C + Ho,m+l + f/irjl)"'^m+l,m , 2 < m < 72 - 2 . (5.9) 

Um^ acts on and leaves 5^m(0) invariant. Our bound bn{C) can then be written 
as 

&n(C) = (^°, (C + Ho,2 + f/i"V'^°)2 • (5.10) 
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The next task is to work out more concretely the operator An,n+i{( + -f^o,n+i) ^ 
^n+i,n appearing in flS.Sp . We use duahty to compute A*, 

{fn,Afn+l)„= / dkfn{ki,...,knyAfn+l{ki,...,kn) 

n „ 

= -2iy/n + 1 / dkdkn+i{sm{27rk(.) +sm{27rkn+i))fn{ki,...,kny 

^/n+l(^l) ■ ■ ■ ,ki — kn+l, • • • , fcji+l) 
» n 

= / dMA;n+i/„+i(/ci, ■ ■ ■ , kn+i) (2iVn+ 1 ( sm{2TT{k(. + k^+i)) 

Jjn+l V ^ 

)* 
(5.11) 

Upon symmetrization, one arrives at 

A*fn{ki,...,kn+i) = 2t^l= y (2sm{2n{kj + ki)) + sm{27ikj) 

+ sin(27rA;^))/„(/ci, . . . ,kj,ke, . . . , kj + ki) . (5.12) 

To compute f/n"'*'^'' we denote the shift ki ki — by Tj. Then 
U^n^^'^Uki, . . . , fc„) = A2(A„,„+i(C + ffo,n+i)"'A„+i,„/„)(A;i, . . . , A;„) 

n „ 

= 4A^ y y / d/c„+i(sin(27r(A;i - kn+i)) + sin(27rA;„+i)) 

x(C + fi„+i(fci, . . .,ke- kn+i, . . . , fcn+i))~Vj(sin(27rA;j) + sin(27rfc<?) 

+2sin(27r(A;j + h)))Tifn{ki, . . . ,kj,ki, . . . , kj + ki) . (5.13) 

Clearly, U^^^'' = Vn"'^^^ + rIt!'~^^\ where Vn"'~^^^ is a multiplication operator and 
rI^~^^^ an integral operator. In the sum in fl5.13p . Vn^~^^^ corresponds to the terms 
i = j and i = n + 1. Hence 

n „ 

Ki"+')(fci,...,A;„) = 2A2V / dfc„+i(sin(27r(fc,-A;„+i)) + sin(27rA;„+i) 

i=i 

+2sin(27r/c^)) {C, + VLn+i{ki, . . . ,ki, - k 

))-'■ (5.14) 

Considering the special cases n = 2, 3 one arrives at the following bounds. 
Theorem 5.1 In the limit — > 0, the following hounds are valid, 

^-12-5/43-3/2^-1/4 < ^^0^ ^ i/)-^^°)° < 2-3/2C-V2 _ (5,15) 
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While the lower bound is not sharp, it establishes that S{j,t) must spread superdif- 
fusively. 



Proof: n = 2: It holds 

1^0x0 ^ /„-0 (/- I zj \-l„-0\0 

, -1 



= [ dA;i(i(2 + cos(27rA;i)))'(C + 2cu(A;i)) 



for C ^ 0. 

n = 3: It holds 



> [C + Ho,2 + X'A,4C + Ho,3y'As,,]-'wy,. (5.17) 



For the quadratic form of A^yl2_3(C + -f^o.s) ^^3,2 one obtains, with /2 G 5^2, 

A'(/2, A3(C + ^0,3)-%,2/2)2 = (4AV3) / dhdk2dh{C + n,ih,k2,h)y' 

(^{2 sm{2n{ki + k2)) + sm{27rki) + sm{27rk2)) f2ih + k2, h) 
+ (2sin(27r(/ci + ks)) + sin(27rA;i) + sin(27rA;3))/2(/ci + fcg, /cg) 
+ (2 sin(27r(A;2 + ks)) + sin(27r/c2) + sin(27rA;3))/2(/c2 + ks, ki)) 

<3(/2,V^2^'72)2, (5.18) 

where Schwarz inequality is used in the last step. By taking limits, the inequality 
holds also for the fiber k = 0. Hence 

{g', (C + H)-'gr > {w', (C + Ho,2 + SV^'Y^wX . (5.19) 

For small ( the dominant part of the integral comes from ki close to and one 
obtains 

{w',{C + Ho,2 + 3Vi)-'wy, 

= [ dfci(C + i2nk,yX^2'/'3^C'^Y' = 2-\X'2'/^3')-'/\-'/^ . (5.20) 

□ 

It is of interest to understand whether Theorem 15.11 could be improved to yield 
the KPZ exponent 1/3 as n ^ 00. Unfortunately, already for n = 4, one would 
need a lower bound on A3,4(C + -^^0,4)^4,3 by a multiplication operator, compare 
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with (IS.lSp . Such a bound does not seem to be available. To make, nevertheless, 
some progress we observe the splitting 

fi-\ = eU^a (5.21) 

m (EH, where \ is a multiplication operator and the remainder ll^-x is an 
integral operator involving either a single particle or a pair of particles. Therefore 
one expects that the limit C will be dominated by the potential v!^}^- In solid 
state physics this type of approximation is called the relaxation time approximation. 
Of course it is uncontrolled, at least at present. 

With this approximation in the expression for o'^^ substitutes vl^\ for U^l^- 
Again for one neglects the integral operator -R|i"^2- Iterating this procedure 

yields the recursion relation for the potentials Vm^ , 

Y^^) = 

Kri (^1) • • • ? ^m) 
m „ 

= 2A^ / dkm+1 (2 sin(27rfc£) + sin(27r(A;£ — fcm+i)) + sin(27rA;m+i)) 

m+1; • • • ; km+l) + ^^+1(^15 ■ ■ ■ : ki km+l, • • • ; ^m+l)) ; 

2<m<n-l. (5.22) 
Denoting 6„(C) within this approximation by dn{C)y one arrives at 

Here w{k) is defined in ( Km . w{0) = 0. 

For small ( the dominant contribution to the integral (15.221) comes from km+i 
close to 0. Let us start with m = n — 1. Then 

V^l\{ku kn-i) = (2A2)91]„_i(A;i, . . . , A;„_i) I dfc„(C + 2o;(A;0)-' 

Jf 

^ 9A2(2C)-'/'^^„-l(fcl, . . . , kn-i) . (5.24) 

The next iteration reads 

^{2X^)9n^.2{ku...,kn-2) f dkn-liC + {9X^)2u{kn-l){20-'^y 

Jt 

- {9\Y\^Cr'^'^n-2{k,, . . . , k^.2) , (5.25) 
and, in general, 

Vi^\ku ...,km) = (9A2)1-"— (2C)-""-'"+^fi„(A;i, ...,km), 

2<m<n-l, (5.26) 
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where the exponents aj are defined recursively through 

Oj+i = |(1 — aj) , «! = . (5.27) 
Substituting fl5.26p in fl5.23p and using w{0) = 1 yields for small ( 

d^{() = 2-1(9A2)-""~i(2C)-"" . (5.28) 
The solution to fl5.27p reads 



1 
3 



a„, = :.(l-(-2)-("-i)), (5.29) 



which for n oo converges to a^o = 1/3. Therefore lim„^oo c^n(C) = c^oo(C) 

cioo(C) = 2-^/^3-2/3(A2c)-i/3^ (530) 

which should be compared with (14.25^ . Remarkably enough, the relaxation time 
approximation yields the KPZ exponent 1/3. The prefactor of (A^C)"^''^ equals 
0.292 in fl4.25p while it is 0.191 in f l5.30p . Thus the relaxation time approximation 
gives a prefactor which is approximately 2/3 off the true value. We take this as 
an indication that one needs more powerful methods to obtain the universal scaling 
form for the two-point function. 



6 Continuum limit 

As remarked in the Introduction, the physically meaningful continuum limit of (11. 5p 
must be such as to preserve the scaling (14.70 of the two-point function. Here we want 
to amplify this point and start on the unit lattice. To avoid the various constants, 
we make the specific choice uq = 1/2, Dq = 1. Let us first consider the linear case, 
Aq = 0. Then on the microscopic scale 

ju.it) = l{u,^iit) - 2u,it) + ^i,_i(t)) + 0(t) - • (6-1) 

For the continuum limit we average over a smooth test function g, g : M. ^ W, 
varying on the scale , 6 <^ 1, and take long times 6~'^t with t = 0{1). Then 




where (/> is a mean zero Gaussian field with covariance 

(0(x, t)0(x', t')) = (27r|t - t'l)"^/^ exp [ - (x - x'f/2\t -t'\]. (6.3) 
Equivalently one can view Uj on the lattice with spacing 6 by defining 

v\x,t) = 5-'/\,/si{S-H) (6.4) 
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with [-J denoting integer part and x G M. Then, integrating against smooth test 
functions, the continuum hmit reads 

lim v\x, t) = (j){x,t) . (6.5) 

6— >0 

By (16. ip satisfies 

t) = 5-^{v\x + 5,t)- 2v\x, t) + v\x- 6, t)) 

+ri/2r^(eL./5j(t)-eL(.-5)/5j(t)). (6.6) 

Comparing with (II. Sp . we observe that the bare coefficients are not rescaled. Taking 
the hmit 5 in (16. 6p yields 

dt(P{x, t) = \dl(t){x, t) + d^i{x, t) (6.7) 

with white noise initial data, in agreement with (16. 3p . (16. 5p . 
Next we include the nonlinearity with Aq = 1. Then 

jUjit) = w,{t) - wj^iit) + - 2M,(t) + M,_i(t)) + 0(t) - e,-i(t) (6.8) 

and the stationary two-point function should scale as 

lim 6S{[6-^x], 6~^/H) = r2/3/KPz(t"'/'x) . (6.9) 

In particular, the correct time scale for the continuum limit is 5~^^'^t, t = 0{1). In 
analogy to (16. 5p one introduces the macroscopic field v^{x,t) by 

v\x,t) = 6-'/\,/si{S-'^'t) . (6.10) 

According (16. 8p . v^{x,t) is governed by the evolution 

t) = 5-^ {w\x, t) - w\x - 5, t)) 

+5^'^5-^\ {v\x + 6,t)- 2v\x, t) + v\x - 6, t)) 
+S'/'S-'/'S~' (^L./.j it) - eL(.-.)/.j (t)) . (6.11) 

Comparing with (ll.Sp we conclude that the nonlinearity is left invariant, thus Ab 
remains fixed, while the linear part is scaled down by ^/6. Of course, it is natural 
to conjecture that, as in the linear case, v^{x,t) has a limit as 5 ^ 0. To identify 
the limit one could also use other, better understood models, like the TASEP. Some 
properties are known [5l[T6], but the full limit of v^{x,t) still has to be identified. 

There is one particular case for which analytical results are available [T71 [T2] . In 
fact, it is the problem Burgers wanted to solve [12]. In (II. 2p one sets = and 
studies the decay of the solution for Gaussian white noise initial data, {u{x, 0)u{x', 0)) 
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= {l/8)S{x — x'). We set Ab = 1. The solution to ( 11.20 is well defined in the limit 
z/b = 0. With this meaning, we set i^b = 0. Then the solution u{x,t) is statistically 
self-similar, in the sense that 

u{x, t) = r^'^u{r^'^x, 1) (6.12) 

in distribution, x u{x, 1) is a stationary Markov process with the generator 

LTfiu) = ^f{u) + du'R{u, u') {fiu') - f{u)) (6.13) 

for functions / : M — » R. The second term corresponds to a Markov jump process. 
Thus u{x, 1) = a + x locally, interrupted by downward jumps from u to u' with rate 
R{u,u'). The rate function is computed explicitly and given by 

R(u,u') = (u-u'y^^I(u-u'), (6.14) 
J{u) 

where I and J are given by their Fourier and Laplace transforms in terms of the 
Airy function Ai, 

jf \ ^ n A exp(Mz) 

^("^ = 2^y_.^^^ 2V3Ai(2-V3.) ' (6.15) 
21{u) = {2nur"' + ^ /_2d.exp(«.)(^!^^p^^ + (2.)V2) . (g.ig) 

Plots of the stationary distribution for Lt, and other quantities, can be found 
in [IT]. 
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